Intense slow beams of bosonic potassium isotopes 
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We report on an experimental realization of a two-dimensional magneto-optical trap (2D-MOT) 
that allows the generation of cold atomic beams of 39 K and 41 K bosonic potassium isotopes. The 
high measured fluxes up to 1.0 x 10 11 atoms/s and low atomic velocities around 33 m/s are well 
suited for a fast and reliable 3D-MOT loading, a basilar feature for new generation experiments on 
Bose- Einstein condensation of dilute atomic samples. We also present a simple multilevel theoretical 
model for the calculation of the light-induced force acting on an atom moving in a MOT. The model 
gives a good agreement between predicted and measured flux and velocity values for our 2D-MOT. 

PACS numbers: 32.80.Pj, 33.80.Ps, 42.50.Vk, 07.77.Gx, 03.75.Be 



I. INTRODUCTION 

In the field of atomic physics, experiments and theo- 
retical studies on degenerate mixtures of cold gases have 
been acquiring a growing interest in the last few years. 
The vast majority of works have featured Fermi-Bose 
mixtures, that have now been realized with several com- 
binations of atoms: 6 Li- 7 Li 1], 6 Li- 23 Na 0, 40 K- 87 Rb 
0, 6 Li- 87 Rb 0. Many experiments focused on the de- 
generate Fermi gas, while the bosonic component has 
only been a tool to reach Fermi degeneracy via sympa- 
thetic cooling. Remarkable achievements include BCS- 
BEC crossover 0, El H> El- fermionic Bloch os- 
cillations 01 1 observation of the 3D Fermi surface ^3 • 
Besides, Fermi-Bose systems display a wealth of inter- 
esting phenomena genuinely related to the presence of 
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FIG. 1: Level diagrams for 39 K and 41 K. In the upper part of 
the diagram is shown the natural abundance of both isotopes. 



two species. A f ew, like interspecies Fano-Feshbach res- 
onances HI!. [13. ITU and boson-induced collapse of the 
Fermi gas [l6|, UJ$ > have already been observed but many 
more have been proposed and still await experimental 
confirmation, as e.g. boson-induced superfluidity (lSf 
and mixed phases in optical lattices 0] . 

On the other hand, Bose-Bose mixtures have remained 
relatively unexplored: after the pioneering work on 41 K- 
S7 Rb |2jJ, to our knowledge only one other double-specie 
condensate has been recently produced, namely 174 Yb- 
176 Yb plf. As for the K-Rb mixture, the recent pre- 
cise determination of the interspecies scattering length 
by Feshbach spectroscopy ^ promises a rich phase- 
diagram for the two species loaded in an optical lattice. 
In view of an experiment on a K-Rb Bose-Bose mixture 
we realized and characterized an intense, slow beam of 
cold bosonic potassium isotopes. 

An intense and reliable source of cold atoms repre- 
sents a favorable starting point for every degenerate 
mixture experiment. Two-dimensional magneto-optical 
traps (2D-MOTs) of rubidium have been investigated by 
several groups [23, 13 an d represent one of the bright- 
est sources of slow atoms. With an appropriate double- 
chamber vacuum system, the 2D-MOT provides high 
atomic flux through a small aperture between the two 
chambers. Therefore, it appears perfectly suited for load- 
ing an ordinary 3D-MOT in ultra-high vacuum (UHV) 
environment, the most widely used pre-cooling stage to- 
wards quantum degeneracy of dilute gases. 

In this work we demonstrate an efficient 2D-MOT for 
bosonic potassium isotopes 39 K and 41 K, starting from a 
natural abundance sample. Previously, only_ a 2D-MOT 
of fermionic 40 K has been implemented [l7|. while for 
bosonic isotopes no bright sources have been realized. 
The peculiarity of bosonic potassium is the tight hyper- 
fine structure of the Di line excited state 4 2 P 3 / 2 |2al26|: 
as shown in Fig. ^ the separation between hyperfine lev- 
els is comparable to the transition linewidth, T = 27r x 6.2 
MHz. This feature represents the main physical differ- 
ence with respect to almost all other alkalis (exception 
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FIG. 2: (Color online) Schematics of the 2D-MOT system. 
Two retroreflected transverse cooling beams, with 3:1 elliptic 
section, cross in the center of a Ti structure. A pair of copper 
wired coils generates the quadrupole magnetic field. Near 
the end of the cell a drilled mirror is placed, tilted 45° with 
respect to the longitudinal axis z, along which a 1.5 mm waist 
push beam is propagating, forcing atoms to travel through the 
hole. A plug beam is placed near the hole in the mirror, and 
a probe beam is used to detect the atomic flux in the UHV 
chamber. 



made for Li), and it deeply affects the sub-Doppler cool- 
ing mechanism [2^|. The aim of the present work is to 
address the question of how efficiently 2D-MOT schemes 
work for bosonic potassium. 

The paper is organized as follows. In Sec. ^ we de- 
scribe the experimental setup and some experimental 
techniques used to generate and monitor the cold atomic 
beam. The results of the characterization of the 2D- 
MOT are reported in Sec. IIIII while a theoretical model 
giving a good description of data is explained in Sec. II VI 
In Sec. we conclude with some brief considerations in 
view of future experiments. 



II. EXPERIMENTAL SETUP 

A. Vacuum system 

The vacuum system consists of two chambers, that we 
call 2D-MOT and UHV chamber. The 2D-MOT appara- 
tus is schematically reported in Fig. [21 Four rectangular 
80x35 mm BK7 anti-reflection coated windows, provid- 
ing optical access to four transverse trapping beams, are 
glued on a metallic rectangular frame, machined from 
a titanium block. The atomic vapor is released in nat- 



ural abundance by dispensers, whose injection current 
controls the pressure in the 2D-MOT chamber. The 
vacuum setup design gives the possibility to connect 
to the UHV chamber another identical 2D-MOT cham- 
ber. In this scheme it is possible to load a Rb-K mixed 
3D-MOT starting from two independently controllable 
atomic beams. 

A flexible, 10 cm long, stainless steel bellow connects 
the 2D-MOT and the UHV chambers. A metallic mir- 
ror, tilted 45° with respect to the longitudinal axis, is 
placed near the end of the 2D-MOT chamber and close 
to the bellow. A 1-mm hole is drilled at the center of the 
mirror to let the atoms exit the 2D-MOT chamber. This 
mirror allows us to send a hollow beam, counterpropa- 
gating with respect to the atoms, so that we can turn to 
a 2D+-MOT configuration [23|. 

The 2D-MOT chamber is pumped by a 20 1/s ion diode 
pump, while in the UHV chamber a 55 1/s ion pump, 
combined with the tiny conductance of the mirror hole, 
generates a pressure ratio between the two chambers up 
to 10~ 4 . Inside the bellow, three graphite tubes of in- 
creasing inner diameter (6 to 10 mm) improve the capa- 
bility of pumping alkali metals. 

The quadrupole magnetic field is generated by a pair 
of 90x40 mm rectangular coils, fixed around the pro- 
truding glass windows. This elongated shape yields the 
required almost 2D axial symmetry to the trapping field. 
The magnetic gradient in the x — y plane is set to 17 
Gauss/cm. The current in each coil is independently 
adjusted to compensate for stray magnetic fields along 
the vertical y direction. The compensation of horizontal 
fields along x is provided by a single large coil, placed far 
from the chamber. In this configuration, the longitudinal 
magnetic gradient along z is a factor of 10 smaller than 
in the x — y plane and it will be hereafter neglected. 



B. Laser system 

The cooling and trapping of alkali atoms requires in 
most cases two laser frequencies which are conventionally 
called repumping and cooling. The force exerted upon 
the atoms is dominated by the cooling light, tuned to a 
closed transition, whereas the repumping light is mainly 
needed to counteract the hyperfine optical pumping in 
the ground state, consequent to off-resonant excitation. 
As such, only a small fraction of the laser power suffices 
for the repumping light. 

The situation is completely different for bosonic potas- 
sium, where the tight hyperfine 4 2 P 3 / 2 level spacings are 
comparable to the natural linewidth T. In this case, the 
cooling transition \F = 2) — ► \F' = 3) is not closed since 
the \F' = 1,2) states are excited with similar rates. A 
fast depletion of the \F = 2) ground state toward the 
\F = 1) ensues, the repumping light needs to be fairly 
intense and the cooling force arises from both frequen- 
cies [2^. Hence, the distinction between cooling and re- 
pumping light makes no sense. However, according to a 
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FIG. 3: (a) Typical acquisition of the fluorescence signal. At 
t=0 ms the plug beam is turned on, and after a delay of ~ 8 
ms the fluorescence starts decaying to zero. (b) Measured 
velocity distribution (dots) of atomic beam obtained by a dis- 
crete derivative of (a) and Gaussian fit to data (solid line). 
The fitted peak velocity is 32 m/s with a FWHM of 4.5 m/s. 



widespread convention, we will also call repumping and 
cooling the transitions indicated respectively by lo\ and 
UJ2 in Fig. 

The u>i light is generated by a commercial high-power 
laser diode with back-emission stabilization grating, pro- 
viding 400 mW at 767.5 nm. The frequency stabiliza- 
tion is achieved by locking this laser using modulation- 
transfer saturation spectroscopy in a reference glass cell 
containing potassium vapor. To generate the uj2 radi- 
ation, a beam of about 40 mW is split from the laser 
output, frequency shifted by an acousto-optic modulator 
(AOM) and then injected into a semiconductor tapered 
amplifier that delivers 900 mW. In this way, only the 
repumping laser needs to be frequency locked and the 
frequency difference between the two sources is precisely 
set by the radio-frequency signal fed into the AOM. Each 
of the two main beams is split into several beams with 
independent intensity and frequency control, achieved by 
means of double passage AOM stages. 

After a fiber mode-cleaning stage, a pair of cylindrical 
lenses gives a 3:1 elliptical shape, with a smaller waist 
of 9.4 mm, to two 2D-MOT transverse beams, which are 
circularly polarized and retroreflected after the chamber. 
The maximum laser power for each transverse beam is 
50 mW and 80 mW for u>\ and L02 respectively. 

Along the z axial direction, we shine an additional o~ + 



polarized beam, called push, with a waist of 1.5 mm. To 
achieve the 2D+-MOT configuration, a fraction of this 
beam is split and directed along the — z direction (slowing 
beam). We recall that the slowing beam has a hollow 
profile since it enters the 2D-MOT after reflection upon 
the drilled mirror. 

Another vertical 2 mW beam, called plug, resonant 
with the cooling transition, intercepts the atomic flux 
near the mirror hole: atoms travel unperturbed to the 
UHV chamber when the plug beam is off and are swept 
away when the plug beam is on j2^ . We use this beam to 
perform time-of-flight (TOF) measurements, as described 
below. 



C. Detection technique 

As in Ref. [2^, we analyze the atomic flux by means 
of TOF fluorescence detection. This is accomplished ex- 
posing the atomic flux to a vertical sheet of light (probe 
beam), 30 cm after the mirror hole, and imaging the emit- 
ted fluorescence into a broad area photodiode. The peak 
intensity of the probe beam exceeds 0.5 W/cm 2 , divided 
between lo\ and u>2, both resonant, in a ratio 1:2. 

Switching on the plug beam we interrupt the atomic 
flux; then, from the analysis of the decaying fluorescence 
signal Sp(t) as a function of time, we obtain the longitu- 
dinal velocity distribution p(v) of atoms. Denoting with 
r the time required for an atom with velocity v along z 
to travel the distance L between the plug beam and the 
probe light sheet, we can write the total flux <I> as: 



f t d 

$ = / p(v)dv, where p(v) = —SjrM, (1 

J V dr 



with t = L/v. Here rj is an experimental coefficient ac- 
counting for the calibration of the photodiode and for the 
collection solid angle. 

Fig. Ola) shows a typical acquisition of the TOF sig- 
nal, while the corresponding velocity distribution p(v), 
as derived from Eq. QJ, is displayed in Fig. &b). To 
obtain the total atomic flux <f>, the peak velocity and the 
velocity spread, we fit the data with a Gaussian function. 



III. EXPERIMENTAL RESULTS 

We report the measured flux and velocity distribution 
of 39 K and 41 K atomic beams as a function of experimen- 
tal parameters. We shifted from the first to the second 
isotope by changing the AOMs frequencies and the laser 
frequency main locking point, preserving the geometri- 
cal alignment of the 2D-MOT beams. The analysis on 
41 K was started only after optimization of frequencies 
and power for the 39 K, the former having lower relative 
abundance (6.7%) than the latter (93.3%). 

The region of vapor pressure spanned in the experi- 
ment goes from 2.9x 10 -8 to 3.9x 10~ 7 mbar, correspond- 
ing to the maximum allowed current in the dispenser. We 
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FIG. 4: Measured 39 K atomic flux normalized to unity as a function of transverse 2D-M0T beams detuning (a), and of their 
intensities (b). The detuning values that maximize the atomic flux are Si = — 3.9F and 82 = — 5.8I\ In (c) and (d) we report 
equivalent results obtained for 41 K atomic beam, for which the optimal detunings are 5i = — 2.5r and 82 = — 3F. No data are 
taken in the patterned regions. 



estimate the background gas pressure to be around 10 -9 
mbar. If not otherwise specified, in the experiment we 
set the K gas total pressure to 7xICP 8 mbar. 

The divergence of the atomic beam is measured by 
imaging with a CCD camera the fluorescence emitted in 
the —x direction. The image profile along z corresponds 
to the Gaussian profile of the short axis of probe light 
sheet, while the vertical (y direction) image profile ex- 
tension is limited by the mirror hole. From the known 
width of the probe beam we calibrate the image mag- 
nification and therefore measure the size of the atomic 
beam in the y direction. Given the distance of the probe 
beam from the mirror hole, we calculate a divergence of 
(34 ± 6) mrad. 

A. Trapping beam parameters 

To characterize the 2D-MOT, we first set the intensity 
of both repumping and cooling 2D-MOT beams to their 
maximum, respectively 50 and 80 mW per beam. We 
then make a scan of both frequencies u)\ and 0J2 , search- 
ing for the values that maximize the fluorescence signal, 
hence the total flux of atoms. 



We only present data on the total flux $ because the 
atomic velocities display no significant variations with 
the 2D-MOT parameters. The measured peak velocities 
span a range between 28 and 35 m/s, while the typical 
distribution spread is 4.5 m/s (FWHM). Experimental 
results for 39 K detunings are plotted in arbitrary units 
normalized to the maximum value in Fig. Ufa). The de- 
tunings are defined throughout as follows: 5\ = lj± — lo\2 
and 82 — t^2~^23, where we denote with wff 1 the atomic 
transition |4 2 Si/ 2 , F) -> |4 2 P 3/2 , F'). Note that for these 
measurements the push beam contains only the repump- 
ing component with a power of 6 mW and frequency 
ujp = u>\2 — 5.2r. A more detailed analysis of the push 
beam features is reported in Sec. lIII 01 The detunings op- 
timizing the atomic flux are 5i = — 3.9T and 62 = — 5.8T, 
which correspond to neither of the two configurations re- 
ported in [23. In that work, a MOT is best loaded with 
both lasers detuned below all hyperfine components, in 
a Raman configuration such as {62 — S±\ = A32 where 
A32 = 2ir x 21 MHz is the hyperfine separation between 
the \F' = 3) and \F' = 2) excited levels. Similar results 
were later confirmed in [3(|. We speculate that, since our 
available laser beam intensities are a fifth to a tenth of 
those in j2^ , we are unable to reach such large values of 
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FIG. 5: Measured 39 K gas pressure p in the 2D-MOT chamber 
reported as a function of dispenser current /. The solid line 
is a fit obtained using Eq. P|. 



detunings. 

Then wc fix the detunings and decrease both the cool- 
ing and repumping beams intensity. Fig. Elk) shows the 
39 K atomic flux as a function of the beams intensity: as 
we see, there are no maxima in the explored range. Thus, 
an increase of the laser power should make the atomic 
beam more intense. 

The analysis on 41 K is performed along the same lines. 
Here, the push beam power is set to 6 mW and its fre- 
quency to Up = il>\2 — 4.5T. Experimental results are 
plotted in Fig. Etc, d). As before, while a peak is clearly 
visible in the frequency dependence of the signal, there 
is none in the intensity-dependent plot. Even for 41 K, 
more laser power would enhance the number of atoms in 
the beam. 

The largest contribution to the error on the measured 
total flux comes from the calibration parameter 77, more 
precisely from the evaluation of the solid angle in which 
fluorescence is collected. We estimate this systematic 
uncertainty to be 20%, while the statistical error is 5%. 



B. Vapor pressure 

In order to investigate the dependence of the atomic 
flux upon the potassium gas pressure, we found it con- 
venient to establish a calibration of the vapor pressure p 
against the dispenser supply current I. This is done by 
probing the atomic vapor using linear absorption spec- 
troscopy. We simultaneously record the Doppler absorp- 
tion profiles of the cooling transition in the 2D-MOT 
chamber and in the reference cell used for laser frequency 
locking. The pressure of the latter is inferred from its 
temperature: 2.2 x 10~ 7 mbar at 43 °C. The calibration 
curve of p as a function of I is well-fit by an exponential 
function: 

p = A + exp(I/I -B), (2) 
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FIG. 6: Measured 39 K atomic flux (a) and peak velocity (b) 
reported as a function of gas pressure in the 2D-MOT cham- 
ber. Over the critical value p = 2.1 x 10~ 7 mbar the atomic 
flux starts to be depleted by collisional effects. The solid line 
is a fit obtained using Eq. iQ. 



with A = (3.8 ± 1) ■ 1(T 8 mbar, I = (0.54 ± 0.07) A and 
B = (26.7 ± 1.8). Both the experimental data and the fit 
for 39 K are reported in Fig. [S] 

The maximum allowed current in the dispenser, 6.4 A, 
sets the highest pressure we reach, i.e. 3.9 x 10~ 7 mbar. 
Below the dispenser ignition point, a residual pressure of 
4 x 10~ 8 mbar is likely due to potassium vapor slowly 
released by the chamber walls. 

In Fig. Eta) we report the measured atomic flux as a 
function of the pressure for 39 K. As we can see, the flux 
increases with pressure until a critical point and then de- 
creases. In the case of 39 K this critical point occurs at 
p = 2.1 x 10~ 7 mbar where a flux of 1.0 x 10 11 atoms/s is 
reached. Such a critical point has been observed also for 
Rb [23I I24I l3l| and ascribed to collisions with the back- 
ground vapor. When the inverse collision rate becomes 
of the same order of the cooling time i c , i.e. the time 
required to reach the mirror hole, the probability for an 
atom to reach the mirror hole drops. Assuming this is the 
case also for 39 K, we fit our data to obtain the collision 
rate. The total flux is: 

/•OC 

$ = / y{t c )dt c (3) 
Jo 

where f(t c ) is the cooling time distribution. 

In the collisionless regime f(t c ) is linearly proportional 
to the density, hence the pressure, in the 2D-MOT cham- 



6 



ber. In presence of collisions with the background vapor, 
f(t c ) is depleted for t c > 7, where 7 = up is the colli- 
sion rate, proportional to the pressure p. Indeed, a single 
collision with an atom at room temperature (v t h = 250 
m/s) is sufficient to remove the cold atoms from the 2D- 
MOT velocity capture range. Then, if we assume that 
the cooling time distribution is a Gaussian 



2vr(T 



exp 



(t c - t )' 
2a 2 



Eq. J2J is modified to 
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(4) 



where Erf(x) = ^= J* exp(— t 2 )dt is the usual definition 
of the error function. 

From the numerical simulations illustrated below, we 
obtain the first and second moment of the cooling time 
distribution: to = 5.1 ms and a = 1.7 ms. We then 
use Eq. (0J to fit the data, with two free parameters: 
an overall scaling factor and the collision rate to pres- 
sure ratio k. The results of the fitting is shown by the 
solid line in Fig.[(Ja). The collision rate given by the fit 
7(s~ 1 ) = (8.6 ± 0.6) x p (10~ 8 mbar) is approximately 
one order of magnitude larger than the value reported in 
[32: 0.3 s" 1 at p = 4 x 10" 9 mbar, i.e. 0.75 s" 1 x p (10~ 8 
mbar) . Our value is closer to the collision rates observed 
with rubidium. A meaningful comparison with Ref. |32| 
requires a more detailed investigation on the role of light- 
assisted collisions, assumed negligible for the 2D-MOT. 

Another noticeable effect upon pressure increase is the 
slowing down of the atomic beam, as shown by the plot 
of the peak velocity in Fig. Efb). This behavior is in 
contrast with the case of the simple rubidium 2D-MOT 
described in [24|, where an increase of the gas pressure 
corresponds to the increase of the peak velocity. There- 
fore we attribute the observed effect to the progressive 
absorption of the push beam along its propagation path 
by the atomic gas since, as reported below, the peak ve- 
locity decreases as the push beam intensity is reduced. 

As mentioned before, when not otherwise specified all 
data reported in this work have been taken at p = 7 x 
10 -8 mbar. 



C. Push beam 

We turn to the experimental investigation of the 
atomic beam behavior when the push beam parameters 
are changed. This is a crucial characterization since we 
observe no atomic flux in absence of the push beam, in ac- 
cordance with the findings in experiments with rubidium 
[HHl] where the 2D-MOT without push beam is not ef- 
ficient in our pressure range. Surprisingly, for potassium 
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FIG. 7: Atomic flux $ (a) and peak velocity (b) measured for 
39 K atomic beam as a function of the push beam power. In 



this analysis we set S p = — 5.2V. 



the push beam works best if one uses only repumping 
light, and therefore we define the push beam detuning 
dp = uj p — W12. We defer the discussion of this point to 
the end of the Section. 

First we report the atom flux and peak velocity, 
Figs.Efa) and[7Jb) respectively, for 39 K at fixed detun- 



ing S p = — 5.2r as a function of power. Above 6 mW, 
the atom number is approximately constant, while the 
peak velocity keeps increasing. Therefore, we use a power 
equal to 6 mW (peak intensity of 170 mW/cm 2 ) in order 
to have the maximum flux with a velocity distribution 
still within the capture range of the 3D-MOT in which 
the atoms will be collected. 

No clear enhancements of the flux stem from the addi- 
tion of a counterpropagating hollow beam along —z direc- 
tion (2D+-MOT configuration). Only a 20% reduction in 
mean velocity was found for a a~ beam polarization. 

Given the hyperfine structure of potassium, one would 
expect to increase the efficiency of the push force using 
both cooling and repumping beams, since this avoids the 
depletion of the ground state population. On the con- 
trary, the atomic beam is deteriorated even by a small 
fraction of cooling light LO2 and it is almost extinguished 
when the 0J2 intensity approaches that of u>\. We at- 
tribute this effect to an increase of both the longitudinal 
velocity and the radial temperature of the atoms, when 
pushed by both frequencies. Instead, with a single fre- 
quency the hyperfine optical pumping confines the action 
of the push force within the 2D-MOT volume and the 
atoms drift freely after the mirror hole. 

One may argue that the push beam should then work 
also with the frequency L02 alone. We experimentally 
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TABLE I: Optimal experimental beams power values (P) and detunings (5) used to generate a cold beam of 39 K and 41 K atoms. 
In the right part of the table the total flux <E> values and the corresponding peak velocities are reported. The pressure working 
point is set to 7 x 10 -8 mbar. The magnetic field gradient is 17 Gauss/cm. 
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P [mW] 


Si [r] 


P [mW] 


82 [T] 


P [mW] 


s P [r] 


[atoms/s] 


[m/s] 


39 R 


50 


-3.9 


80 


-5.8 


6 


-5.2 


6.2 x 10 10 


35 


41 R 


47 


-2.5 


85 


-3.0 


6 


-4.5 


5.2 x 10 9 


33 



find that, at equal power, 0J2 light is far less efficient 
than cji. We speculate that this is due to lo 2 being 
more prone to hyperfine optical pumping for two rea- 
sons: (i) in the 2D-MOT there is more cooling than re- 
pumping light, (ii) UJ2 detuned a few linewidths to the 
red of the \F = 2) — > \F' = 3) transition is near to the 
\F = 2} — ► \F' = 1, 2) resonances causing hyperfine op- 
tical pumping (the same effect is less important for the 
uji light, since the transitions \F = 1) — > \F' = 2, 1) 
are closer and there is no hyperfine optical pumping for 
\F = 1) — > \F' = 0)). Another possible way to avoid hy- 
perfine optical pumping is to use a blue detuned cooling 
light |33|. even if this would increase the pushing effi- 
ciency mainly on fast atomic classes and hence confer to 
the atomic flux a higher peak velocity. 

As Fig. E[a) shows, decreasing the intensity of the laser 
beam the atomic flux decreases as well. In the pure 
2D-MOT configuration, when no push beam is shone on 
atoms, no flux is detected, independently of both power 
and frequency of transverse beams and quadrupole mag- 
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FIG. 8: Measured atomic flux $ and peak velocity of 39 K (a, 
b) and 41 K (c, d) as a function of the push beam frequency 
detuning S p . The push beam power is set to 6 mW. 



netic field gradients. 



D. Summary 

In Tab. [J are summarized all the experimental param- 
eters which maximize the total atomic flux $ of the two 
isotopes and the values of the corresponding fluxes. For 
p = 7xl0~ 8 mbar high atomic fluxes are achieved, con- 
taining 6.2 x 10 10 atoms/s and 5.2 x 10 9 atoms/s for 
39 K and 41 K respectively. The pressure value and their 
small peak velocities, respectively 35 and 33 m/s, are 
perfectly adequate to an efficient subsequent loading in 
a UHV environment 3D-MOT. By increasing p up to 
2.1 x 10~ 7 mbar a flux of 1.0 x 10 11 atoms/s of 39 K is then 
observed. Nearly the same gain is expected for 41 K flux 
as well at the same value of pressure p. 

The high intensity of atomic fluxes is an essential point 
in order to reach a degenerate regime for K isotopes fol- 
lowing a double MOT cooling and trapping scheme. At 
optimal detunings, with a 15 times lower light intensity, 
we achieve a 3D-MOT loading rate higher than the one 
reported in j^. Thus, we conclude that, to generate 
a cold atomic beam, a 2D-MOT shows a much higher 
efficiency when compared to an ordinary 3D-MOT. 



IV. A SIMPLIFIED MULTILEVEL 
THEORETICAL MODEL 

To have a quantitative description of our observations 
we use an extension of the theoretical model discussed in 
[3l| . This model shortcuts the integration of the optical 
Bloch equation by assuming a heuristic expression of the 
total force exerted by the different beams, all having the 
same frequency, on a two-level atom: 



h 



(5) 



I s T 2 + 4(5, - k, • v) 2 



where i,j denote the beams and I s = irhcT/(3A 3 ) is 
the two-level saturation intensity, equal to 1.8 mW/cm 2 
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for potassium D2 line. Authors of Ref. [3lJ employ this 
model to analyze a 3D-MOT of Rb with the addition of 
a push beam. The extension of this treatment to bosonic 
potassium, because of the narrow upper level structure, 
requires to take into account all the allowed hypcrfinc 
transitions. For this purpose, we introduce the further 
assumption that forces arising from different transitions 
add independently. This consistently disregards coher- 
ences among the two 4Si/ 2 hyperfme states, which how- 
ever play no role in the Doppler cooling mechanism. In 
principle, we should consider even the Zceman structure 
of the hyperfme levels; in practice, to reduce the number 
of transitions contributing to the total force, we calculate 
the detunings and the line strength in a manner to aver- 
age out the Zeeman sublevels. For each laser beam i and 
each transition \F) — > \F'), we define the average detun- 
ing Af F as the center-of-mass of all the Zeeman com- 
ponents weighted by the squared Clebsch-Gordan coeffi- 
cients \c(F', m'\ 1, o"i, F, m)\ 2 , where <Ji denotes the beam 
polarization (all beams are circularly polarized with the 
MOT required helicity). This detuning depends linearly 
on the displacement from the 2D-MOT axis via the mag- 
netic field gradient. We then define the strength of each 
hyperfme transition: 



Gf F ' = 



E m , m '\<F'^m';l,a u F,m)\ 2 



We incorporate hyperfme optical pumping by break- 
ing the force into two parts, due to the cooling and the 
repumping light, weighted by the relative populations in 
the F = 2 and F = 1 ground levels. Therefore, the 
expression (0 of the total force is generalized as follows: 



Pi 



-fi 



P2 



Pi +P2 Pl+ Pi 



(0) 
(7) 



with 



Si,F 



Ij,F \ -« n FF 



r 2 +4(^-k 4 -v-Af F ')' 



where Sf — s j- Here, Ii tl ,I% ia are respectively the 
repumping and cooling intensity of the i-th beam, and 
61,62 the corresponding detunings as defined earlier. The 
populations p\ , P2 are taken as the equilibrium values of 
the rate equations for the six hyperfme components. 



A. Numerical Simulations 

A numerical integration of the classical equations of 
motion yields the phase-space trajectories. We consider 
only the atoms that, at t = 0, lie on the boundary surface 





(d) 



2 4 6 S 

Push beam power (mW) 



■12 -10 -8 -6 -4 -2 

Push beam detuning (fi /r) 



FIG. 9: Results of simulations for 39 K : total flux and peak 
velocity as a function of the push beam power (a,b) and de- 
tuning (c,d). 



S of the 2D-MOT volume, approximated by a rectangu- 
lar box with sizes equal to the beam diameters. In sam- 
pling the velocity-space, the Boltzmann factor is nearly 
unity for all velocities lying within the 2D-MOT capture 
range, well below the 250 m/s thermal velocity. From the 
integration we extract: (i) the fraction p of the atoms ex- 
iting the mirror hole to the atoms entering the 2D-MOT 
volume, (ii) the longitudinal velocity distribution of the 
atomic beam, and (Hi) the distribution of the cooling 
time <p(t c ), as defined earlier. 

To obtain the total flux, we only need to multiply p by 
the total number of atoms entering the cooling volume 
per second, at pressure p and room temperature T: 



<j> = px5 



P 



(8) 



In our simulation multiple scattering of light and intra- 
beam atomic collisions are neglected. The collisions with 
background gas, occurring at rate 7 = 60 s _1 , are ac- 
counted for by weighting each trajectory with a factor 
exp(— -ft c ) to deplete the tail of atoms with long cooling 
times t c . We also select only those atoms flying in a cone 
34 mrad wide around the longitudinal axis z. 

Setting the experimental parameters as in Table|H with 
a quadrupole magnetic field gradient of 17 Gauss/cm, for 
39 K we obtain a total flux <I> sim = 8.7x 10 10 atoms/s, with 
a peak velocity of 30.2 m/s, in good agreement with the 
experimental values reported in Table [I] 

Repeating the simulation for different values of the 
push beam power and detuning we find the curves plotted 
in Fig. |5| The agreement with the experimental findings 
is satisfactory for the peak velocities. As for the total 
flux, the model shows saturation in push power at lower 
values than in the experiment and an optimum detuning 
of — 6r close to the observed value of — 5.2r. The calcu- 
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lated dependence of atomic flux on push beam power and 
detuning reproduces only qualitatively the experimental 
curves of Figs. 6(a) and 7(a). We believe that coherent 
effects ignored in our simulations are likely responsible 
for the discrepancies with data. A more exact analysis 
based on the integration of optical Bloch equations is 
needed to address these issues but is beyond the scope of 
our simplified model. 

V. CONCLUSIONS AND OUTLOOK 

The main conclusion of this work is that, similar to ru- 
bidium, potassium atoms can be coaxed into an intense 
and cold beam, still preserving the UHV constraints. Re- 
cently, we have used two separate 2D-MOT beams to load 
a double-specie K/Rb MOT in the UHV chamber as a 
preliminary step towards quantum degenerate Bose-Bose 
mixtures. The combination of potassium and rubidium 
offers a wealth of possibilities to tune the interspecies 
interactions via Feshbach resonances [22^. Broad Fano- 
Feshbach resonances are predicted for 39 K- 87 Rb and 41 K- 



87 Rb in the absolute ground state at 318 G and 515 G, 
respectively. For 39 K, sympathetic cooling with Rb can 
be a route to single condensates with tunable interac- 
tions, where resonances are expected around 40 G for two 
atoms in the magnetically trapped state \F — 1, m = — 1) 
[sil and around 20 G for the absolute ground state [35l |. 

An exciting perspective is to tune interactions in the 
presence of an optical lattice, where a rich diagram of 
insulating and superfluid phases is predicted when the 
interspecies scattering length is varied 36]. Also, as 
suggested in [$3, entering the Mott-insulator regime for 
both species paves the way to Feshbach-associated long- 
lived bosonic heteronuclear molecules, since three-body 
losses can be suppressed by tailoring the Mott phase to 
have one atom per species in each lattice site. 

This work had financial support from Ente Cassa di 
Risparmio in Fircnze and from Istituto Nazionale di 
Fisica Nucleare through project SQUAT. We warmly 
thank P. Clade for his help, G. Modugno for many fruitful 
suggestions, C. Fort and all the members of the Quantum 
Degenerate Gas group at LENS for discussions. 
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